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Abstract 



Extending to 4 degrees of freedom a symplectomorphism used in attitude dynam- 
ics it is shown in a direct way the connection between the 4-D isotropic harmonic 
oscillator and the 3-D Kepler systems. This approach made transparent that only 
when we refer to rectilinear solutions, the bilinear relation defining the KS trans- 
formation is needed. 



1 Introduction 



More than 40 years after the KS transformation pQ, the epitome of which might be the 
work of Cordani [2], any proposal of a new insight on the connection between the 4-D 
isotropic oscillator and Kepler systems might be taken as the opposite. Nevertheless, 
the way in which many authors still deal with the relation between these integrable 
systems, has pushed finally to take the risk and write this letter. Our claim is that, 
in order to present the connection of these systems we do not need to go by concepts like 
embeddings, weak extended canonical transformation formalism j3] , or the more recently 
used geometric algebraic approaches [U [5j [6]. We came to understand what is explained 
here studying the 4-D isotropic oscillator reductions [TJ. More precisely, it is shown in a 
direct way that the Kepler- Coulomb flow appears as a part of the 4-D isotropic oscillator 
flow. 

The core of our approach goes back to a basic concept in 3-D dynamics 'the instanta- 
neous plane of motion ' and the associated nodal-polar variables, already used in planetary 
theories by Hill [8j and much later in satellite theory (see Deprit [9]). Today some au- 
thors refer to them as Whittaker transformation HDl.ITT]. Nevertheless, for reasons which 
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are unknown to this author, this chart still continues to be ignored, versus the spherical 
variables, not only in vast field of mechanics but even in the special area of Hamiltonian 
astrodynamics. 

Inspired in the 3-D Whittaker transformation, what we propose here is to consider 
part of the Andoyer angles (see f[T2"j) ) as the ones defining a 'plane of motion' joint with 
a function of the distance, which lead to the Kepler system as part of the 4-D oscillator. 
This author is convinced that the use of 'polar-nodal variables' in 4-D context will be of 
great benefit also in quantum mechanics, molecular dynamics, etc. 

The 4-D isotropic oscillator is an integrable dynamical system defined by the para- 
metric Hamiltonian function 

^ = ^£(Q* W), (i) 

where a; is a parameter. There is a large and uninterrupted literature about this system 
going from Jauch and Hill [12] to Waldvogel [13], including Moser jTJ] and Iwai [15] , due 
to the fact that it is one of the rare few examples of maximally superintegrable systems 
(see Fasso [To]). 



2 Starting with Projective Euler variables 

We consider the transformation: VEp : (p,(fi,9,ip) — > (<7i, (fe ?3> (Z4), dubbed as Projective 
Euler variables, given by 

qi = F(p) sin - cos 9s = F(p) cos - sin (2) 

<?2 = F(p) sin - sin — ^— , q 4 = F(p) cos - cos — — , 

with (p,(j),8,tp) G R + x [0, 27r) x (0,7t) x (— |,fj- When F(p) = 1, the transformation 
defines Euler parameters as functions of Euler angles. We choose here F(p) = J~p. This 
transformation is well known in the literature (see, for instance, [T7J [3j HHJ [191 120] ) . 
The canonical extension associated to the transformation ([2]) is obtained as a Mathieu 
transformation, which satisfies J2Qidqi — P dp + Q d(f> + Q d9 + ^ dip. The relations among 
the momenta are given by 

P = * 2 (giQi + ?2<?2 + qsQs + qiQi), 

^ i-i qi 

Q = (giQi + q2Q2) {qj + qf> - (q^Qs + q^Q^iql + ?I) 



^(ql + qlM + ql) 
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$= ^(-^Qi + q^ + qiQ-s-q-sQi), (3) 

Note that the factor y/p has a long history. Indeed, as Bartsch [6] put it: 'For the one- 
dimensional Kepler motion, it was already found by Euler [2T] that the introduction of 
a square-root coordinate u = y/x and a fictitious time r defined by dt = xdr reduces 
the Kepler equation of motion to d 2 u/dr 2 + 2Eu = i.e. the equation of motion of a 
one- dimensional harmonic oscillator'. This reemerges in Heggie and Hut [22] (pp. 145) 
who seem unaware of this work of Euler. 

Then, excluding the invariant manifolds Ai\ = {{q,Q)\qi = q-i = 0} and M. 2 = 
{(<Z>Q)|<73 = <?4 = 0} where Levi-Civita transformation already shows the connection of 
the 2-D Kepler and isotropic oscillators, the Hamiltonian ([Q) in the new variables may be 
written as 



P U , o„o2 , 2 4,2 , $ 2 + ^ 2 - 2$^ COS^' 

sin 2 ^ 




(4) 



i.e. variables and ip are cyclic, with $ and \l/ as the corresponding first integrals. 
Note that a dash is used instead of the variable to stress the fact that this coordinate is 
ignorable. In other words the differential system reduces to 

dp _dU d9 _ dU dP dU dO &R 

J?~~dP' d^'de' ~d^~~~dp' ~d7~~~de 

and two quadratures 

(f) = J(dH/d$)dT and ij) = J(dH/d^)dr. (5) 

In order to reach the main result of this letter, let us remember first two basic features of 
our system, presented in the following Propositions, based on the reordering of the terms 
defining Eq. (j3J) and the suitable chosen regularizing function. 

Proposition 1.- The 4-D isotropic oscillator may be considered as a generalized Kepler 
system. It reduces to a 3-D Kepler system and a quadrature when we restrict to either 
$ = or \l/ = manifolds. 

Proof - Let us consider Poincare technique with a time regularization r — > s given by 
dT=(A P y l ds. (6) 
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Then, the flow is defined by the Hamiltonian Iti = j-(H — h) where h is a fix value of the 
Hamiltonian H for chosen initial conditions, and the flow is defined now on the manifold 
7^ = 0. In a slight different form, it is more convenient to write 



h 



2 V P 
in the manifold /C = - 



1 [ p i © 2 $ 2 + ^ 2 -2$^ cos#\ 

2 p 2 sin 2 9 J 4p 



(7) 
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Therefore, denoting part of function (J7j): 

H-K = 77 \ P 2 + -7f + . 2 J — , 

2 \ p 2 p 2 sin 9 J p 
then, the Hamiltonian function (J7J) may be also written as 
$2 _ 2 $# C os# 



JC = U K + 



2p 2 sin 2 



(9) 

) in spherical 
(10) 



Note that the function Hk is the Hamiltonian of the Kepler system (see 
coordinates 

x = p sin 9 cos 0, y = p sin 9 sin <fi, z = p cos # 
(where 9 is the colatitude) and their momenta (P, 0, $), if we choose 
7 = h/4, 



where 7 is the fundamental constant of the Kepler- Coulomb system. Observe that we 
might have taken the term \l/ 2 /(p 2 sin 2 9) instead of $ 2 /(p 2 sin 2 #) for the definition of 
T-Lk- 



The differential system defined by (Q is given by 



dp 


dt 


dU K 


ds 


dP ~ 


dP ' 


d9 


dt 


dU K 


ds 


ae " 


dQ ' 


dP 


dt 


dU K 


ds 


dp 


dp 
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&H K 
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d9 



*(* - 2$ cos#) 
p 3 sin 2 9 

*(^cos#-$(l + cos 2 fl)) 
p 2 sin 3 9 
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and the two quadratures (J5]). Then, if we restrict to the manifold \1/ = and we identify 
the variable s with the physical time t, the flow of the oscillator given by Eqs. (ITTi) joint 
with the quadrature <fi = J(d7iK/d^) ds, reduces to a Keplerian flow in the manifold 
%k = — f • The last equation is the quadrature of if) 



f dK. , f $ cos 6 

to be computed after the Kepler system is integrated. q.e.d. 

This shows that the 4-D oscillator may be seen as a generalized Kepler system. In fact, 
in the literature this Hamiltonian relates to Hartmann and other ring-shaped potentials 
but we will not enter that issue here. 

Note that, considering the inverse of the Projective Euler transformation fl2]), we may 
see the transformation from spherical to Cartesian ffTO]) as a projection R A — > R 3 . Explic- 
itly, inverting we have p = J2qj an d 



■ a 2A n g3 + gj ~ 9l ~ g2 

smO = — , cos = — . 



P 



P 



sm. 



sin-0 



<?1<?3 + <?2<?4 , qm - g2<?3 

cos (b — — 



A 

qm - 92^4 
A 



cosip 



A 

qm + <?2<?3 



where A = J (ql + q^iqi + ql), we obtain immediately 
x = 2(gig 4 - q 2 q 3 ), 

y = 2(gig 3 + 92^4), 
2,2 2 2 

z = q 3 + q i -q 1 - q 2 , 



in other words, the KS-transformation. 



3 Switching to Projective Andoyer variables 



The Andoyer variables (introduced by Serret [25], and also referred as Serret- Andoyer [26] 
or Andoyer-Deprit [27J) are a well known symplectomorphism in dynamical astronomy 
|28[ |29| [30] and recently introduced in other fields such as attitude and control [31] . The 
reader ought to be aware that in the literature authors use different letters for them (see 
[28| I29| 132]). In what follows the (A, fi, u, A, M, N) notation is adopted. 
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Assuming the vector ($, 0, \l/) different from zero, i.e. excluding the invariant man- 
ifold of rectilinear solutions treated before, the canonical transformation from Andoyer 
(A, fx, v, A, M, N) to Euler (0, 9, ip, $, 6, is given by 



cos e cos a — sin e sin a cos p — cos 9 = 0, 

cos 9 cos a + sin 9 sin a cos — v) — cos e = 0, 

cos 9 cos e + sin # sin e cos (0 — A) — cos a = 0, 

$ = A, (12) 
* = A^, 



Q= I 2 _A? + N*-2NAc O b0 



sin 2 9 

where cose = A/M and cos a = N/M. In fact, the transformation requires two charts. 

When we add to them the variables (p, P) we obtain what we call the 'Projective 
Andoyer' transformation: 

f p, (j), 9, ip \ ( p, A, //, z/ 



P, $, 8, * ^ i F, A, M, JV ) ^ 



which, versus Projective Euler transformation, this one is not a canonical extension. 
Note that completing the expressions of the momenta as functions of {q,Q), after some 
computations we obtain 

1 



M=- v /||g|| 2 ||g|| 2 -(g-g) 2 . (14) 

which shows that when Projective Andoyer variables are not defined: M = 0, the motion 
is rectilinear. 

Proposition 2.- In Projective Andoyer variables, the system defined by the 4-D isotropic 
oscillator, properly regularized, is separable in two subsystems, one linear in the angle p 
and a quadrature for the variable p. 

Proof- Expressing our Hamiltonian (J3J) in the 'Projective Andoyer' (PA) variables it 
results 

up : n _ n2 : 2M 2 



n u = H(p, -, -, -, P, -,M,-) = -f + 2pP 2 + 



2 p 

We see that the three angle variables (X,p, v) are cyclic. In fact, not only the momenta 
(A, N) but also the variables (A, v) are integrals. 
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Then, fixing a value of the Hamiltonian h and making a change of independent variable 
dr = (p/4)ds the Hamiltonian takes the form 

r p iu m up 2 p 2 P 2 hp 1 2 

in the manifold /C = 0. From this Hamiltonian we obtain immediately the conclusions. 
Details are not needed for our purposes. q.e.d. 

Now we present the main result of this letter: 

Theorem.- In Projective Andoyer variables, the system defined by (2P, properly regular- 
ized, includes the Keplerian system for any value of the integral N. 

Proof- After a change of independent variable according to Poincare technique, /C = 
g(p)(H u — h) considering now ([H]), i.e. g(p) = l/(4p), the Hamiltonian of the 4-D isotropic 
oscillator is given by 

in the manifold /C = —u/8. Now, what remains is to connect this Hamiltonian with the 
Kepler system in 3-D. The use of the polar-nodal canonical transformation (see Deprit 
[TO] ) extended to four dimensions 

(p, A, p, u, P, A, M, N) -> (x, y, z, u, X, Y, Z, N) (16) 

is one way to do it. Indeed, denoting lZ(v, a) a rotation matrix of angle a around 
the vector v, the transformation ([TBI) is defined considering three rotations related to 
the three direct orthogonal reference frames (ei,e2,es), (^1,^2,^3) and (61,62,63), where 
£1 = K(e 3 , A)ei, 63 = I)e 3 and 61 = 1Z(b 3 , p)£ u by 

(x,y,z) T = n(e 3 ,\)n(£ 1 ,I)n(b 3 ,p)(p,0,0) T 

(x,y,z) t = n(e 3 ,\)n(£ u i)n(b 3 ,p)(p,—,o) T , 

p 

with cos / = A/M, and where T stands for transpose of a vector. Explicitly, we have 

x = p(cos p cos A — sin p sin A cos I) , 
y = p(cos p sin A + sin p cos A cos /) , 
z = p sin p sin/, 
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and similarly for the momenta. It is easy to verify that ( |T6l) is a canonical transformation 
of Mathieu type: Xdx + Ydy + Zdz = Pdp + AdX + Mdp, but not a canonical extension. 
Perhaps for this reason is not so well known. Then, our Hamiltonian (ITS"]) expressed the 
variables (x, y, z, z/, X, Y, Z, N) takes the form 

U = H{x, y, z, -, X, Y, Z,-) = \ \\X\\ 2 - (17) 

which is the Hamiltonian of the 3-D Keplerian system where, with some abuse of notation, 
we have written ) and X = (X, Y, Z). Moreover, (u, N) are integrals that take 

any value. q.e.d. 

Let give details of the inverse process: Let consider a reference frame {ei,e2,e^). 
We define the vectors u and n such that x = pu, where p = \\x\\ and = 1, and 
M = x x X = Mn, with M = ||M||, ||n|| = 1. Then, we take the momenta P, A as 
P = x ■ and A = xY — yX. Moreover, we may write x n = (sm\)£, with 

\\£\\ = 1, and cos A = e\ ■ £, join with cos p = £ ■ u, and smp = (n x £) • u, which 
ends the inversion of the polar nodal transformation. Then from Eqs. f lT2|) we obtain the 
projective Euler canonical variables. Finally, from Eqs. ([2]) we obtain qi. It rests to have 
the expression for Qi (i — 1,4). From the canonical extension Eqs. ([3]) we know they are 
linear in them. Thus, we obtain explicitly Qi by inverting the matrix associated with the 
transformation, but it is not necessary to be given here. 

We see that, comparing with KS and in contrast with it, no constraint is needed for the 
integrals (z/, N). Then, what gives KS transformation such a special place in the oscillator 
Kepler connection? The fact that Projective Andoyer transformation is not defined when 
M = 0, which according to Eq. ( |T1"]) corresponds to rectilinear trajectories: q\ \ Q. 

Note that, although the Projective Andoyer variables do not yet define a set of action 
angle variables, they are an intermediary step in that direction. Indeed, from (JTSJ) we 
readely obain a Delaunay set of variables for 4-D oscillators, a symplectomorphism well 
suited for perturbation theories. We just need to make use of the Delaunay transformation 
(p,p,P,M) ->• (8,g,D,G) (see [23], [10]). It reduces Eq. ([15]), i.e. (Op), to only one action 
U = - 1 2 /(2D 2 ) (for details [33]), which reflects the fact that our systems are maximally 
superintegrable [IB] . 

(Qi, <?2, Q3, Q4\ 

Wl, Q2, Qs, QJ 



( x, y, z, v \ 
V X, Y, Z, N J 

This diagram, based on the group structure of canonical transformations, suggests to 



Proje ctive Euler / p, (J), 6, ijj \ 

~^ v p, $, e, * ) 

PA 

Polar Nodal / P, A, //, V \ 

^~ \ P, A, M, N ) 
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carry out explicit compositions of some of them, just to obtain more insight, as well as 
possible relations with other known transformations, or new ones that might be defined 
likewise. Further relations among the previous transformations, the option F(p) = p 
included, as well as the connection with other integrable and perturbed systems is in 
progress [33] . 

Thanks to Dr. Egea for his help in an early version of this letter and to Dr. Lara for his 
comments. Partial support came from projects MTM2006-06961 and MTM2009-10767 of 
the Ministry of Technology and Science of Spain, and a grant from Fundacion Seneca, of 
the Autonomous Region of Murcia. 
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